System coordinator and modular architecture for open-loop and closed-loop control of diabetes

ABSTRACT

A structure, method, and computer program product for a diabetes control system provides, but is not limited thereto, the following: open-loop or closed-loop control of diabetes that adapts to individual physiologic characteristics and to the behavioral profile of each person. An exemplary aspect to this adaptation is biosystem (patient or subject) observation and modular control. Consequently, established is the fundamental architecture and the principal components for a modular system, which may include algorithmic observers of patients&#39; behavior and metabolic state, as well as interacting control modules responsible for basal rate, insulin boluses, and hypoglycemia prevention.

CROSS-REFERENCES TO RELATED APPLICATIONS

The present application claims priority from U.S. application Ser. No. 13/322,943, filed Nov. 29, 2011, entitled “System Coordinator and Modular Architecture for Open-Loop and Closed-Loop Control of Diabetes;” of which the disclosure is hereby incorporated by reference herein in its entirety.

FIELD OF THE INVENTION

Some aspects of some embodiments of this invention are in the field of medical methods, systems, and computer program products related to managing the treatment of diabetic subjects, more particularly to glycemic analysis and control.

BACKGROUND OF THE INVENTION

People with diabetes face a life-long optimization problem: to maintain strict glycemic control without increasing their risk for hypoglycemia [13,58,59]. The engineering challenge related to this problem is to design algorithms using automated insulin delivery to exert optimal closed-loop control of glucose fluctuations. Since the early studies of continuous external glucose regulation (e.g., BioStator, [10]), two primary approaches have emerged: The use of classic proportional-Integral-derivative (PID) algorithms, and modem methods based on models of the human metabolism. The first studies using subcutaneous insulin delivery and continuous glucose monitoring (CGM) employed PID control [57,60]. Recently, model predictive control (MPC), received considerable attention [20, 21, 44, 50] due to its many clinical and engineering advantages.

BRIEF SUMMARY OF THE INVENTION

MPC is typically based on a model of the human metabolic system. Fortunately, the modeling of glucose-insulin interaction is one of the most advanced applications of mathematics to medicine. Beginning with the now classic Minimal Model of Glucose Kinetics (MMGK) co-authored by Dr. Claudio Cobelli who leads the Italian team of this project [2], a number of elaborate models have been developed [16,21]. These models can be classified in three broad classes: (i) models to measure parameters that are not accessible by direct lab tests, such as MMGK assessing insulin sensitivity; (ii) models to simulate that enable in silico pre-clinical trials, and (iii) models to control used to empower algorithms such as MPC.

An aspect of an embodiment of the present invention provides the progress towards advisory open-loop control or automated closed-loop control that will be greatly accelerated by a structured modular approach to building control components. Specifically, an aspect of an embodiment of the present invention provides a system of control modules responsible for basal rate, pre-meal and correction insulin boluses, and hypoglycemia prevention. These modules will be informed by biosystem observers providing information about the patients' glycemic state. A modular approach to closed-loop control development would have a number of advantages that include, but are not limited to:

-   -   Incremental testing of modules in parallel or consecutive         studies;     -   Incremental FDA approval and industrial deployment of system         features;     -   User flexibility—each system observer or control module could be         used separately, or within an integrated control system,         depending on patients' or physicians' choice;

An aspect of an embodiment provides an external open-loop or closed-loop control that shall have separate interacting components responsible for prevention of hypoglycemia, postprandial insulin correction boluses, basal rate control, and administration of pre-meal boluses. These control modules receive information from biosystem observers that are responsible for tracking glucose fluctuations and the amount of active insulin at any point in time. This dual control-observer architecture is dictated by the natural separation of the computational elements of a closed-loop control system into algorithms observing the person and algorithms actuating control. Central role in this architecture is played by the system coordinator—an algorithmic module that is responsible for controlling the integration and the interactions of the modular system.

Pertaining to the feasibility of each of the proposed observers and control modules, as well as the feasibility of using CGM technology and subcutaneous insulin delivery for automated closed-loop control the following may be referenced:

-   -   Use and accuracy of CGM; algorithmic processing of CGM data: [3,         5, 24, 26, 29, 34, 40];     -   In silico pre-clinical trials: [15, 16, 28, 31];     -   Glucose variability observer and Risk Analysis: [36, 37, 38, 39,         42, 43, 46];     -   Insulin observer, subcutaneous insulin transport, sensitivity,         and action: [0, 2, 7, 12, 47, 54];     -   Control module 1: prediction and prevention of hypoglycemia: [4,         26, 35, 53, 56, 61];     -   Control modules 2, 3, and 4: correction boluses and closed-loop         control: [8, 11, 27, 44, 45, 51].

An embodiment of the present invention defines a modular architecture that can accommodate a variety of system observers and control modules that can be assembled into a system for open-loop advisory mode control or automated closed-loop control of diabetes.

An aspect of an embodiment of the present invention, a structure, method, and computer program product for a diabetes control system provides, but is not limited thereto, the following: open-loop or closed-loop control of diabetes that adapts to individual physiologic characteristics and to the behavioral profile of each person. An exemplary aspect to this adaptation is biosystem (patient) observation and modular control. Consequently, an aspect of an embodiment of the present invention establishes the fundamental architecture and the principal components for a modular system, which includes algorithmic observers of patients' behavior and metabolic state, as well as interacting control modules responsible for basal rate, insulin boluses, and hypoglycemia prevention. An exemplary role in this architecture is played by the system coordinator—such as an algorithmic module that may be responsible for controlling the integration and the interactions of the modular system.

An aspect of an embodiment of the present invention provides a structure, method, and computer program product for a diabetes control system provides, but is not limited thereto, the following: open-loop or closed-loop control of diabetes that adapts to individual physiologic characteristics and to the behavioral profile of each person. An exemplary aspect to this adaptation is biosystem (patient) observation and modular control. Consequently, an aspect of an embodiment of the present invention establishes the fundamental architecture and the principal components for a modular system, which may include algorithmic observers of patients' behavior and metabolic state, as well as interacting control modules responsible for basal rate, insulin boluses, and hypoglycemia prevention.

An aspect of an embodiment of the present invention provides a structure for a diabetes control system. The structure may comprise: modules for processing and storing data; conduits between modules; and signals produced in the event that certain modules are not inserted within the structure.

An aspect of an embodiment of the present invention provides a computer program product comprising a computer useable medium having a computer program logic for enabling at least one processor in a computer system for a diabetes control system. The computer program logic may be configured to include: modules for processing and storing data; conduit means between modules; and producing signals in the event that certain modules are not inserted within the system.

An aspect of an embodiment of the present invention provides a method for enabling a diabetes control system. The method may comprise: providing modules for processing and storing data; providing conduits or the like between modules; and producing signals in the event that certain modules are not inserted within the system.

These and other objects, along with advantages and features of various aspects of embodiments of the invention disclosed herein, will be made more apparent from the description, drawings and claims that follow.

BRIEF DESCRIPTION OF THE DRAWINGS

The accompanying drawings, which are incorporated into and form a part of the instant specification, illustrate several aspects and embodiments of the present invention and, together with the description herein, serve to explain the principles of the invention. The drawings are provided only for the purpose of illustrating select embodiments of the invention and are not to be construed as limiting the invention.

FIG. 1 schematically provides an exemplary embodiment of the structure of the modular architecture for control of diabetes.

FIG. 2 schematically provides an exemplary decoupled basal and bolus control supervised by hypoglycemia prevention.

FIG. 3 schematically provides a first exemplary embodiment of the general architecture for control of diabetes.

FIG. 4 schematically provides a detailed embodiment of the general architecture for control of diabetes.

FIG. 5 schematically provides an exemplary embodiment of the general architecture for control of diabetes.

DETAILED DESCRIPTION OF THE INVENTION Definitions

For the purposes of an embodiment of the present invention:

-   -   Open-loop advisory-mode control is defined as a system that uses         CGM data and information from insulin pump to provide real-time         advice about treatment adjustments to patients with diabetes;     -   Closed-loop control system is defined as a system using data         from CGM devices and information from insulin pump to         automatically control insulin delivery by the insulin pump.

Overview of Modular Architecture of Open-Loop and Closed-Loop Control of Diabetes:

As shown in FIG. 1 , a modular open- or closed-loop control system includes observers of patients' metabolic state, daily profiles and behavior, which provide information to the central system coordinator, which in turn directs the actions and the interactions of an array of control modules. Specifically:

A System Coordinator will coordinate the distribution of input signals to control modules (routing) and most importantly allocate different segments of diabetes management to different controllers by restricting the input to these controllers (see example in discussion below). Finally, the system coordinator will ensure the direct feeding of external inputs to controllers if the observers are inactive.

Observers will receive frequent information about metabolic measurements (such as continuous glucose or insulin), metabolic disturbances (such as meals or exercise), and metabolic treatments (such as insulin or glucagon injections). Based on these inputs the observers will

-   -   construct and update an internal representation of the metabolic         state of the patient and transmit this state to the control         modules;     -   keep an internal representation of the behavioral pattern of the         patient, such as daily meal and exercise profiles; and     -   assess risks for undesirable events such as hypoglycemia or         glucose variability.

Safety Supervision Module will receive information from the observers and from the control modules (below) and will decide whether there is an increase of the risk for upcoming hypoglycemia or prolonged hyperglycemia. If risk increase in encountered, the module will reduce or discontinue the suggested insulin infusion.

Three Control Modules will be responsible for insulin administration. The Control Modules will receive instructions from the System Coordinator and will supply their output to the Safety Supervision Module for evaluation. The Control Modules are:

-   -   Control Module 1 will calculate and suggest the basal insulin         delivery;     -   Control Module 2 will calculate and suggest compensation of the         basal delivery (up or down) in case of non meal related         deviations, such as drops or rise due to exercise, or residual         dawn phenomenon not covered by module 1;     -   Control Module 3 will calculate and suggest meal insulin         boluses, potentially including pre-meal priming boluses.

While the Safety Supervision Module and the Control Modules are subjects of independent invention disclosures or have been developed elsewhere [44,51], the subject of an embodiment the present invention is defining the general architecture of an open-loop or a closed-loop control system and the interactions between the components (modules) of that system as presented in FIGS. 1 and 2 . A feature of an embodiment of the present invention is the design of the System Coordinator, responsible for the seamless integration of observers, control modules, and safety supervision.

In the context of engineering design of an embodiment, it may be important to underscore that:

-   -   Each system observer or control module can be used separately,         or within an integrated open- or closed-loop control system,         depending on patients' or physicians' choice. This modular         approach will allow the incremental testing of system features         in parallel or consecutive studies;     -   The operation of the system in open-loop advisory mode will be         conceptually similar to closed-loop control mode, but will         differ in terms of implementation: open-loop will provide         real-time information to the patient, while closed-loop will         control directly the insulin pump.

An exemplary idea behind the introduction of a System Coordinator is its ability to decouple different control functions, and to coordinate the control module action with separate modules responsible for different aspect of diabetes management, such as meals, exercise, basal pattern, and hypo/hyperglycemia avoidance. In other words, separate interacting algorithms will suggest optimal pre-meal bolus control (e.g. a stochastic algorithm) and will exercise basal rate control or administer post-meal correction boluses (e.g. deterministic algorithms). This “separation of duties” corresponds to the stochastic nature of meals and behavior, and the deterministic nature of basal and postprandial physiology, and also has deep mathematical reasoning motivated by the experience gained in our recent clinical trials of closed-loop control in Type 1 diabetes. These trials showed excellent overnight regulation but rather slow (as compared to open loop) breakfast regulation. Tuning of the algorithm aggressiveness alone was not sufficient to achieve both goals. It was therefore necessary to introduce a strategy that handles differently night and breakfast regulation. With this in mind, we introduce the System Coordinator, which allows each Control Module to operate within a certain BG range. Specifically, after Module 3 treats a meal, the effect of this bolus and the concurring meal is projected 1-2 h ahead and will be subtracting in real time from the trajectory from the incoming CGM data. In other words, the System Coordinator will “correct” the CGM track sent to Module 2 by the projected action of Module 3. This modus operandi is illustrated in FIG. 2 : the observed glucose trace is presented by a gray line, but Module 2 will only “see” the black trace, which is the difference between real-time CGM and the action of Module 3. This will result in a rather simple interaction between Modules 1, 2 and 3: Module 1 sets the reference treatment for the day, Module 3 will operate in a regime with updates every several hours, its action will be continuously projected, and the result will be supplied to Module 2, which will operate in frequent increments (e.g. 15 minutes) Further, all control modules will be supervised by the Safety Supervision module, which will warn the person for upcoming hypo/hyperglycemia and will suggest reduction in insulin delivery or correction boluses in open-loop advisory mode, or will directly reduce or discontinue the insulin pump infusion rate in closed-loop control mode.

FIG. 3 presents a flow chart of the interactions between the system components identified in FIG. 1 , while FIG. 4 presents a detailed view of the nodes and the conduits of the modular architecture, featuring the interactions of the System Coordinator with the other components of an open-loop advisory system or a closed-loop control system. The only difference between open-loop advisory operation and closed-loop control operation is in the delivery of information from the System Coordinator: in open-loop the information about insulin rate and potential for hypoglycemia is presented to the patient; in closed-loop control mode insulin delivery commands are send directly to the insulin pump. The following paragraphs describe the function of each of the modules of the architecture, focusing largely on the interfaces between modules and the requirements/definition of each signal. This description assumes that the system updates are computed in discrete time t, corresponding to given sampling interval. We use the terms “time,” “discrete time,” and “stage” interchangeably.

The Data Module serves to process raw data from external sources, scan the data for integrity, and produce four real-time signals that are then used by the remaining modules of the modular system. As shown in FIG. 4 , the external data sources are (1) continuous glucose monitoring (CGM) data, possibly from multiple sensors, (2) data from the patient's insulin pump (i.e. insulin actually delivered to the patient), (3) exercise information (e.g. heart rate information and other indicators of physical activity), and (4) meal information (e.g. acknowledgements of meals as they arrive). Only the CGM and the insulin pump data are mandatory; exercise and meal information is optional and will be used only in certain embodiments of this invention. The outputs of the data module are:

-   -   g_(c)(t)=a single, processed glucose sample at stage t,         reflecting the effect of all processing of the glucose sensors         used as input     -   J(t)=the most recent (relative to t) actual insulin pump command         (mU/min) It should be noted that all insulin injections are         expressed as rates of infusion, including boluses.     -   EX(t)=exercise process data at time t (EX is an optional time         series that is specific to the embodiment of the modular         architecture)     -   M(t)=meal data at stage t (M is an optional time series that is         specific to the embodiment of the modular architecture)

The Long Term Observer Module serves to produce statistical profiles of the patient's meal and exercise behavior, β and η respectively, and a circadian profile the patient's daily insulin utilization γ, which are used by other control modules to set priors on parameter estimates and/or constraints on insulin action. The inputs to the Long Term Observer Module are all of the Data Module outputs shown in FIG. 4 : g_(c)(t), J(t), EX(t), and M(t). The outputs of the module are:

-   -   β(τ)=a daily profile of the patient's meal behavior, where the         notation (τ) expresses the fact that the profile is used as a         lookup table descriptive of the patients behavior as a function         of time of day     -   η(τ)=a daily profile of the patient's exercise behavior, where         again the notation (τ) expresses the fact that the profile is         used as a lookup table descriptive of the patients behavior as a         function of time of day     -   γ(τ)=a daily (circadian) profile of the patient's utilization of         insulin, where again the notation (τ) expresses the fact that         the profile is used as a lookup table descriptive of the         patients behavior as a function of time of day         All outputs of the Long Term Observer Module are optional the         entire Long-Term Observer would only exist in certain specific         embodiments of the modular architecture. Some embodiments of the         modular architecture will not include an active Long Term         Observer; in these embodiments the outputs of the model take         default values, β⁰(τ) η⁰(τ), and γ⁰(τ), which are not utilized         by the remaining modules of the system.

In one embodiment of the Long Term Observer Module, the meal behavioral profile β defines a probabilistic description of the patients eating behavior within given “meal regimes” throughout the day. Specifically, for each meal regime β would comprise the conditional probability of a meal arriving within the next sampling interval (within the meal's window of opportunity), given that the meal has not yet arrived:

$\begin{matrix} {p_{k} = \frac{f_{k}}{1 - F + {{\sum}_{i = 1}^{\overset{\_}{k}}f_{i}}}} & (1) \end{matrix}$

where f_(k) is the frequency with which the corresponding meal arrives in the k-th sampling interval of the current meal regime, 1-F is the probability with which the meal will not arrive within its window of opportunity, and k is the latest possible stage in which the meal could arrive. In this scheme it is historical information about meals M(t) that allows for the “observation” of the meal behavioral profile β(τ), expressed as conditional probabilities p_(k). Note that if the meal time is known in advance within the j-th sampling interval, then p_(k) will be zero for all k not equal to j, and p_(j)=1. This embodiment of the meal behavioral profile would allow for the administration of insulin in anticipation of meals, without compromising patient safety.

The Short Term Observer Module of FIGS. 3 and 4 serves to compute real-time estimates of key metabolic states for the patient, including possibly plasma insulin and glucose concentration, which are used by Control Modules 1-3 to perform various tasks. The Safety Supervisor, for example, may use metabolic state estimates to assess/predict the risk of out-of-range excursions in blood glucose (i.e. episodes of hypoglycemia or hyperglycemia).

The biometric values that comprise the “key” metabolic states are specific to the embodiment of the modular architecture. However, at least one of the states will be blood glucose concentration. We use {tilde over (x)}(t) to denote the vector of true metabolic state values, reserving the first element of the vector {tilde over (x)}₁(t) to be defined as blood glucose concentration. It is important to note that the true values of the metabolic states are unknown in general because they are not the same as the input signals g_(c)(t), J(t), EX(t), and M(t). (For example, while plasma glucose and insulin cannot be measured in real time, values for these states can be estimated from g_(c)(t), J(t), EX(t), and M(t).) We use {circumflex over ({tilde over (x)})}(t) to denote the corresponding vector of state estimates.

Thus, the inputs of the Short Term Observer Module are the same as the outputs of Data Module shown in FIG. 3 : g_(c)(t), J(t), EX(t), and M(t). The output of the Short Term Observer Module is:

-   -   {circumflex over ({tilde over (x)})}(t)=a vector of estimates of         the key metabolic states of the patient, including plasma         glucose and plasma insulin.

In embodiments of the modular architecture where the Short Term Observer is not active, the default state estimate is simply a pass-through of the glucose sample:

{circumflex over ({tilde over (x)})}(t)=g _(c)(t).  (2)

Even though the composition of the vectors {tilde over (x)}(t) and {circumflex over ({tilde over (x)})}(t) is specific to the embodiment, we can describe the basic framework for how state estimation is performed. We assume that the evolution of the state vector is described by a discrete-time, nonlinear dynamic model, generally expressed as:

{tilde over (x)}(t+1)={tilde over (F)}({tilde over (x)}(t),J(t),ω_(m)(t),ω_(c)(t)),  (3)

where ω_(m)(t) and ω_(c)(t) are disturbance processes representing meals and exercise/physical activity. The nonlinear model could be a discrete-time realization of the oral glucose meal model [16], minimal models derived from the meal model, the Hovorka [21] or the Sorensen model [55], or another yet to be determined mathematical model of glucose-insulin-exercise dynamics. The state estimate {circumflex over ({tilde over (x)})}(t) is derived through a process of filtering (state “observations”) based on the processed input data g_(c)(t), J(t), EX(t), and M(t). The filtering process could be a direct application of Kalman filtering, extended Kalman filtering, or another yet to be determined statistical procedure that incorporates state observation as a key internal process. The state estimation process is driven in part by a model for the relationship between these signals and the underlying state vector, expressed as:

CGM(t)={tilde over (G)} _(CGM)({tilde over (x)}(t))+ν_(CGM)(t)  (4)

M(t)={tilde over (G)} _(m)(ω_(m)(t))+ν_(m)(t)  (5)

EX(t)={tilde over (G)} _(c)(ω_(c)(t))+ν_(c)(t),  (6)

where (1) {tilde over (G)}_(CGM), {tilde over (G)}_(m), and {tilde over (G)}_(c) describe the functional relationship between the subject's actual metabolic state {tilde over (x)}(t) and CGM(t), the actual meal disturbance ω_(m)(t) and the meal signal y_(m)(t), and the actual exercise/physical activity disturbance ω_(e)(t) and the exercise data signal y_(e)(t), respectively and (2) ν_(CGM)(t), ν_(m)(t), and ν_(e)(t) are sensor noise processes.

The System Coordinator shown in FIG. 4 plays a central role in the modular architecture, serving (1) to act as a “router” for signals from and to the various Control Modules and (2) to decompose the estimate blood glucose concentration into parts attributable to (i) meals and the response to meals and (ii) all other disturbances (exercise and physical activity). The System Coordinator receives a specific set of signals from the various Control Modules, listed in the following paragraphs.

From the Data Module, the System Coordinator receives g_(c)(t), J(t), EX(t), and M(t). From the Long Term Observer Module, the System Coordinator receives the behavioral profiles β(τ), η(τ), and γ(τ), and, from the Short Term Observer Module, the System Coordinator receives the metabolic state vector estimate {circumflex over ({tilde over (x)})}(t). In addition, from the Safety Supervision Module, the System Coordinator receives:

-   -   u_(m,actual)(t)=insulin allowed at time t by the Safety         Supervisor as a part of the response to meals     -   u_(r,actual)(t)=insulin allowed at time t by the Safety         Supervisor as a part of the response to non-meal disturbances

From the Meal Control Module (Control Module 3), the Coordinator Receives

-   -   g_(m)(t)=estimated glucose excursion due to meals (and the         response to meals)

Many of the inputs to the System Coordinator are “passed through” to other control modules, as shown in FIG. 3 .

In addition to acting as a “router” for signals between various modules, the System Coordinator play a central role in attributing glucose excursions to either (1) meals and the response to meals or (2) other metabolic disturbances. Specifically, the System Coordinator serves to evaluate

-   -   g_(r)(t)=an estimate of blood glucose concentration offset by         the estimated contribution of meals (and the response to meals),         computed as

g _(r)(t)={circumflex over ({tilde over (x)})}(t)−g _(m)(t)  (7)

The signal g_(r)(t) is the only output of the System Coordinator that is not a pass-through signal from other modules. (Note that since output g_(m)(t) of the Controller defaults to zero if the Meal Control Module is not active, then g_(r)(t) is exactly the estimated value of blood glucose {circumflex over ({tilde over (x)})}(t).)

The Daily Profile Control Module (Control Module 1) of FIG. 4 serves to compute a reference insulin signal u_(ref)(t), which is analogous to the patient's daily insulin profile in conventional CSII therapy. The inputs to the Daily Profile Control Module are the outputs of the system observers: (1) the metabolic state vector estimate {circumflex over ({tilde over (x)})}(t) (produced by the Short Term Observer Module) and (2) the behavioral profiles β(τ), η(τ), and γ(τ) (produced by the Long Term Observer Module). The output of Daily Profile Control Module is:

-   -   u_(ref)(t)=reference insulin infused at time t

The computation of the reference insulin signal is specific to the embodiment of the modular architecture. In an embodiment where the Daily Profile Control Module is inactive, the default output is simply the patient's daily basal insulin profile.

The Meal Control Module (Control Module 3) of FIG. 4 serves to manage BG excursions in anticipation of and in response to meals in a quasi open-loop fashion. It does this by (1) maintaining an open-loop internal estimate of the contribution of meals to the patient's overall BG and (2) using the patient meal behavioral profile to compute recommended meal insulin u_(m)(t) both prior to the meal arriving and after the meal. The inputs to the Meal Controller are

-   -   M(t)=the meal data process (from the Data Module)     -   u_(m,actual)(t)=insulin allowed at time t by the Safety         Supervisor as a part of the response to meals (produced by         Safety Supervision Module)     -   β(i)=the meal behavioral profile (produced by the Long Term         Observer Module) The outputs of the Meal Controller (originating         within the Meal Controller) are     -   u_(m)(t)=the recommended insulin infusion at time t for         accommodating meals     -   g_(m)(t)=estimated glucose excursion due to meals (and the         response to meals)

The computation of u_(m)(t) and g_(m)(t) are both specific to the embodiment of the modular architecture. In an embodiment where the meal control module is inactive, the recommended insulin infusion u_(m)(t) defaults to zero, as does g_(m)(t) the estimated glucose excursion due to meals, leaving the Compensation Control to Target Module (if active) to reject the meal disturbance.

For embodiments in which the Meal Control Module is active, estimation of the glucose excursion due to meals would be computed from an open loop dynamic system model, reflective of the dynamic interactions between glucose and insulin. Whatever dynamic system model is used can generally be described as a nonlinear discrete-time system, with state vector 4(t), whose evolution is dictated by:

ξ(t+1)=F _(m)(ξ(t),u _(m,actual)(t),M(t)).  (8)

where the state space equations defined by the operator F_(m) are such that, if u_(m,actual)(t) is fixed at zero and M(t) is also fixed at zero, then ξ(t) converges to zero asymptotically. The estimated excursion in blood glucose can be derived from the model as

g _(m)(t)=G _(m)(ξ(t)).  (9)

The same dynamic system model could be used in conjunction with the meal behavioral profile β(τ) to compute an optimal open loop response to meals, perhaps even in anticipate of meal arrival.

In one embodiment of the modular architecture with an active Meal Control Module, the behavioral profile β(τ) could express as the conditional probability p_(k) of the meal arriving at stage k of the current meal regime (given that it has not already arrived). Such a meal behavioral profile would allow, the administration of insulin in anticipation of meals, without compromising patient safety. In this embodiment, the suggested meal control signal would be computed as a function of the control relevant statistics y_(m)(t) and future conditional probabilities p_(t), p_(t+1), K, p _(k) as

u _(m)(t)=C _(m)(ξ(t),M(t),p _(t) ,p _(t+1) ,Kp _(k) )  (10)

where C_(m) denotes the mathematical transformation of input data to optimal anticipatory and reactive insulin injections.

The Compensation Control to Target Module of FIG. 4 (Control Module 2) serves to manage in-range excursions of blood glucose that are not already accounted for by the Meal Control Module (Control Module 3). In serving in this capacity, this module maintains an internal estimate of the patient's “residual” metabolic state based on (1) the meal-offset glucose signal g_(r)(t) produced by the System Coordinator, past residual insulin injections u_(r,actual)(t) allowed by the Safety Supervisor, the exercise data process EX(t), and the exercise behavioral and insulin utilization profiles η(τ) and γ(τ).

The inputs to the Compensation Controller are

-   -   g_(r)(t)=an estimate of glucose excursions away from g_(ref)(t)         that are not due to meals (computed by the System Coordinator)     -   u_(r,actual)(t)=insulin allowed at time t by the Safety         Supervisor (computed by the Safety Supervisor)     -   EX(t)=the exercise data process from the Data Module)     -   η=exercise behavioral profile (produced by the Long Term         Observer Module)     -   γ=circadian insulin profile (produced by the Long Term Observer         Module) The output of the Compensation Controller is     -   u_(r)(t)=the recommended “residual” insulin infusion at time t         needed to compensate for non-meal disturbances

The computation of u_(r)(t) is specific to the embodiment of the modular architecture. In an embodiment where the Compensation Control to Target Module is inactive, the default is u_(r)(t)=0. In this case the Safety Supervision Module (in conjunction with the Daily Profile Control Module and Meal Control Module) would take full responsibility for keeping the patient within an acceptable range of BG values.

In an embodiment of the modular architecture with an active Compensation Control to Target Module, the recommended residual insulin could be computed using Proportional-Integral-Derivative (PID) control, as in [57], using e(t)=g_(r)(t)−g_(ref)(t) as an error signal, where g_(ref)(t) is a possibly time-varying target.

In other embodiments of the modular architecture (in which the Compensation Control to Target Module is Active), the recommended residual insulin signal could be computed from a dynamic model of the interactions of g_(r)(t) and residual insulin, which in general terms could be described as:

x(t+1)=F _(r)(x(t),u _(r,actual)(t),ω_(e)(t)),  (11)

where the state space equations defined by the operator F_(r) are such that, if u_(r,actual)(t) is fixed at zero and the exercise disturbance is also fixed at zero (i.e. no physical activity), then x(t) converges to a given target value g_(ref) asymptotically. A vector of state estimates {circumflex over (x)}(t) can be derived (through a process of state observation) using measurement models

g _(r)(t)=G _(r)(x(t))  (12)

EX(t)=G _(e)(ω_(e)(t))+ν_(e)(t),  (13)

The recommended residual control signal u_(r)(t) is computed based on an estimate {circumflex over (x)}(t) of the state vector x(t) that is computed from the residual glucose signal g_(r)(t). In one embodiment of the modular architecture, recommended residual control signal would be computed using closed-loop model predictive control techniques, as in [44] based on the estimate {circumflex over (x)}(t). In other embodiments, Compensation Control could be achieved via LQG [51], LMPC, or any other closed loop control methodology. In yet other embodiments, compensation control could be achieve via either positive or negative residual “boluses”, whose timing and extent are computed based on {circumflex over (x)}(t).

The Safety Supervision Module of FIG. 4 monitors and, if necessary, modifies the suggested insulin injection signals produced by the other control modules, in an effort to (1) avoid hypoglycemia and (2) if the other modules are unable to do so, mitigate sustained hyperglycemia. In addition, the safety supervisor is responsible for computing the final “approved” next insulin injection J(t+1) used as the main output of the modular architecture, attributing relevant parts of this signal to an “actual” meal control signal u_(m,actual)(t) and an actual residual control signal u_(r,actual)(t).

The inputs to the Safety Supervisor are

-   -   {circumflex over ({tilde over (x)})}(t)=the vector of estimated         metabolic states for the patient, the first component of which,         {circumflex over ({tilde over (x)})}₁(t), is an estimate of the         patients blood glucose concentration (Short Term Observer         Module)     -   M(t)=the meal data process (from the Data Module)     -   EX(t)=the exercise data process (form the Data Module)     -   u_(ref)(t)=reference insulin infused at time t (produced by the         Daily Profile Control Module (Control Module 1))     -   u_(m)(t)=the recommended insulin infusion at time t for         accommodating meals (produced by the Meal Control Module         (Control Module 3))     -   u_(r)(t)=the recommended residual insulin infusion at time t for         accommodating non-meal disturbances (produced by the         Compensation Control to Target Module (Control Module 2))         The outputs of the Safety Supervisor (all generated by the         Safety Supervisor) are     -   u_(m,actual)(t)=insulin allowed at time t by the Safety         Supervisor as a part of the response to meals     -   u_(r,actual)(t)=insulin allowed at time t by the Safety         Supervisor as a part of the response to non-meal disturbances     -   J(t+1)=total insulin infusion allowed by the Safety Supervisor         at time t The computation of J(t+1) and the decomposition of         J(t+1) into actual meal and residual components is specific to         the embodiment of the modular architecture. It is a requirement,         however, that all of approved insulin must be account for as         reference, meal, and residual insulin:

J(t+1)=u _(ref)(t)+u _(m,actual)(t)+u _(r,actual)(t)  (14)

In one embodiment of the modular architecture, the computation of J(t+1) would be based on the prediction of future metabolic states {circumflex over ({tilde over (x)})}(t+τ|t), given the current estimate {circumflex over ({tilde over (x)})}(t) of the metabolic state vector, the exercise data process EX(t), and the suggested insulin amounts u_(ref)(t), u_(m)(t), and u_(r)(t) held fixed at their current values:

{circumflex over ({tilde over (x)})}(t+τ|t)=Φ_(r)({circumflex over ({tilde over (x)})}(t),EX(t),u _(ref)(t),u _(m)(t),u _(r)(t)).  (15)

The predicted future metabolic states {circumflex over ({tilde over (x)})}(t+τ|t) would be factored into an assessment of the risk of hypoglycemia {circumflex over (R)}(t+τ|t) using the risk symmetrization procedure of [36,43]. In this embodiment, the total insulin control signal would be computed from the suggested values as

$\begin{matrix} {{J\left( {t + 1} \right)} = \frac{{u_{ref}(t)} + {u_{m}(t)} + {u_{r}(t)}}{1 + {k \cdot {\hat{R}\left( {{t + \tau}❘t} \right)}}}} & (16) \end{matrix}$

where k is an aggressiveness factor for the attenuating function of the module, which would be determined from patient characteristics, such as body weight, total daily insulin, carb ratio, etc. The difference between J(t) and u_(ref)(t)+u_(m)(t)+u_(r)(t) would be attributed to u_(m,actual)(t)+u_(r,actual)(t) according to a proportionality constant, α in [0,1], as follows:

u _(m,actual)(t)=α·(J(t+1)−u _(ref)(t))  (17)

u _(r,actual)(t)=(1−α)·(J(t+1)−u _(ref)(t))  (18)

Note that if the risk {circumflex over (R)}(t+τ|t) of hypoglycemia is zero (which would be the case when predicted BG is greater than 112.5 mg/dl), then J(t+1)=u_(ref)(t)+u_(m)(t)+u_(r)(t), u_(m,actual)(t)=u_(m)(t), and u_(r,actual)(t)=u_(r)(t).

An analogous computation could be used to compute a total approved insulin J(t) to mitigate hyperglycemia, attributing the difference to u_(m,actual)(t) and u_(r,actual)(t).

Turning to FIG. 5 , FIG. 5 is a functional block diagram for a computer system 500 for implementation of an exemplary embodiment or portion of an embodiment of present invention. For example, a method or system of an embodiment of the present invention may be implemented using hardware, software or a combination thereof and may be implemented in one or more computer systems or other processing systems, such as personal digit assistants (PDAs) equipped with adequate memory and processing capabilities. In an example embodiment, the invention was implemented in software running on a general purpose computer 500 as illustrated in FIG. 5 . The computer system 500 may includes one or more processors, such as processor 504. The Processor 504 is connected to a communication infrastructure 506 (e.g., a communications bus, cross-over bar, or network). The computer system 500 may include a display interface 502 that forwards graphics, text, and/or other data from the communication infrastructure 506 (or from a frame buffer not shown) for display on the display unit 530. Display unit 530 may be digital and/or analog.

The computer system 500 may also include a main memory 508, preferably random access memory (RAM), and may also include a secondary memory 510. The secondary memory 510 may include, for example, a hard disk drive 512 and/or a removable storage drive 514, representing a floppy disk drive, a magnetic tape drive, an optical disk drive, a flash memory, etc. The removable storage drive 514 reads from and/or writes to a removable storage unit 518 in a well known manner. Removable storage unit 518, represents a floppy disk, magnetic tape, optical disk, etc. which is read by and written to by removable storage drive 514. As will be appreciated, the removable storage unit 518 includes a computer usable storage medium having stored therein computer software and/or data.

In alternative embodiments, secondary memory 510 may include other means for allowing computer programs or other instructions to be loaded into computer system 500. Such means may include, for example, a removable storage unit 522 and an interface 520. Examples of such removable storage units/interfaces include a program cartridge and cartridge interface (such as that found in video game devices), a removable memory chip (such as a ROM, PROM, EPROM or EEPROM) and associated socket, and other removable storage units 522 and interfaces 520 which allow software and data to be transferred from the removable storage unit 522 to computer system 500.

The computer system 500 may also include a communications interface 524. Communications interface 124 allows software and data to be transferred between computer system 500 and external devices. Examples of communications interface 524 may include a modem, a network interface (such as an Ethernet card), a communications port (e.g., serial or parallel, etc.), a PCMCIA slot and card, a modem, etc. Software and data transferred via communications interface 524 are in the form of signals 528 which may be electronic, electromagnetic, optical or other signals capable of being received by communications interface 524. Signals 528 are provided to communications interface 524 via a communications path (i.e., channel) 526. Channel 526 (or any other communication means or channel disclosed herein) carries signals 528 and may be implemented using wire or cable, fiber optics, blue tooth, a phone line, a cellular phone link, an RF link, an infrared link, wireless link or connection and other communications channels.

In this document, the terms “computer program medium” and “computer usable medium” are used to generally refer to media or medium such as various software, firmware, disks, drives, removable storage drive 514, a hard disk installed in hard disk drive 512, and signals 528. These computer program products (“computer program medium” and “computer usable medium”) are means for providing software to computer system 500. The computer program product may comprise a computer useable medium having computer program logic thereon. The invention includes such computer program products. The “computer program product” and “computer useable medium” may be any computer readable medium having computer logic thereon.

Computer programs (also called computer control logic or computer program logic) are may be stored in main memory 508 and/or secondary memory 510. Computer programs may also be received via communications interface 524. Such computer programs, when executed, enable computer system 500 to perform the features of the present invention as discussed herein. In particular, the computer programs, when executed, enable processor 504 to perform the functions of the present invention. Accordingly, such computer programs represent controllers of computer system 500.

In an embodiment where the invention is implemented using software, the software may be stored in a computer program product and loaded into computer system 500 using removable storage drive 514, hard drive 512 or communications interface 524. The control logic (software or computer program logic), when executed by the processor 504, causes the processor 504 to perform the functions of the invention as described herein.

In another embodiment, the invention is implemented primarily in hardware using, for example, hardware components such as application specific integrated circuits (ASICs). Implementation of the hardware state machine to perform the functions described herein will be apparent to persons skilled in the relevant art(s).

In yet another embodiment, the invention is implemented using a combination of both hardware and software.

In an example software embodiment of the invention, the methods described above may be implemented in SPSS control language or C++ programming language, but could be implemented in other various programs, computer simulation and computer-aided design, computer simulation environment, MATLAB, or any other software platform or program, windows interface or operating system (or other operating system) or other programs known or available to those skilled in the art.

An embodiment of the present invention provides for, but not limited thereto, the development of an open-loop advisory system assisting patients with diabetes in the control of their blood glucose level, or a closed-loop control system (known as artificial pancreas) controlling automatically blood glucose. Accordingly, an aspect may provide, but not limited thereto, a complex project involving the development and the implementation of multiple mathematical algorithms, methods and engineering solutions.

Some aspects of an embodiment of the present invention diabetes control system, method and computer program product provides, but not limited thereto, the following: open-loop or closed-loop control must adapt to individual physiologic characteristics and to fast changing environmental factors; and keys to this adaptation are biosystem (patient) observation and modular control. Consequently, an aspect of an embodiment of the present invention diabetes control system, method and computer program product establishes the foundation for a modular system comprised of algorithmic observers of patients' behavior and metabolic state, and control modules responsible for insulin delivery and hypoglycemia prevention.

A component of such a modular system is the System Coordinator, which integrates the action of observers, control, and safety modules in order to:

-   -   In open-loop advisory mode, advise the patient about changes in         basal rate, need for insulin boluses, or need to discontinue         insulin delivery and take a preventive action due to increased         risk for hypoglycemia.     -   In closed-loop control mode, command directly the insulin pump         and issue warning for impending hypoglycemia to the patient if         insulin pump shut-off would be insufficient to eliminate that         risk.

Moreover, an embodiment of the present invention includes a modular approach enabled by a System Coordinator that will permit incremental testing and deployment of system features—observers and control modules—which will structure and facilitate system development.

Existing open- and closed-loop control algorithms do not include modular architecture—typically a single control module is implemented and charged with the function of delivering insulin regardless of the causes of glucose fluctuation. A modular approach enabled by a System Coordinator of an embodiment of the present invention allows the distribution of control and safety functions among specialized control modules that are then integrated by the System Coordinator.

A modular architecture of open-loop advisory system or closed-loop control system and method has many advantages, including but not limited to: the possibility of incremental development, testing, and deployment of control modules, and the possibility of using existing modules for incorporation in the system. Centralizing the system integration and coordination functions into a System Coordinator alleviates the separate control modules from the need of working in regimes that are unsuitable for their design. For example, a deterministic MPC is poorly suited to account for the stochastic nature of carbohydrate intake during meals, while a stochastic impulse control is poorly suited for maintaining a steady basal rate. Employing a System Coordinator permits different types of specialized algorithms to function together, each within the realm of its optimal performance.

The Modular Architecture and the System Coordinator included in an embodiment of the invention is suitable for implementation in open-loop advisory systems or closed-loop control systems for diabetes. These systems would typically use a continuous glucose monitor and an insulin pump, linked by the System Coordinator to optimize glucose control in diabetes. Other sources of information (e.g. heart rate monitoring that enables the recognition of exercise) can be included in the system as well, as long as interface with the System Coordinator is established.

It should be appreciated that as discussed herein, a subject may be a human or applicable in principle to animals. It should be appreciated that an animal may be a variety of any applicable type, including, but not limited thereto, mammal, veterinarian animal, livestock animal or pet type animal, etc. As an example, the animal may be a laboratory animal specifically selected to have certain characteristics similar to human (e.g. rat, dog, pig, monkey), etc. It should be appreciated that the subject may be any applicable human patient, for example.

Examples

Practice of an aspect of an embodiment (or embodiments) of the invention will be still more fully understood from the following examples, which are presented herein for illustration only and should not be construed as limiting the invention in any way.

An aspect of an embodiment of the present invention provides a structure, system, or method for a diabetes control system. An embodiment of the structure or related method may comprise: modules for processing and storing data; conduits (or the like) between modules; and signals produced in the event that certain modules are not inserted within the structure. The structure may comprise different implementations of modules that can be inserted and/or interchanged. The modules for processing and storing data may be configured to include one or more of the following: one or more data acquisition modules; one or more observer modules; one or more routing modules; one or more control modules; and one or more safety modules. The structure may further comprise an insulin injector for injecting insulin based on the output of the one or more safety modules. One or more of the data acquisition modules may receive one or more of the following types of information: continuous glucose monitoring data; insulin pump data; exercise data; and meal data. The exercise data may include heart rate information, motion sensor information, and other indicators of physical activity. The meal data may include acknowledgements of meals as they arrive.

The one or more data acquisition modules may be configured to output one or more of the following types of information: a single, processed glucose sample at a specific time, or a history of glucose samples up to a specific time, or a statistic computed from glucose samples up to a specific time; a most recent actual insulin pump command, or a history of recent insulin pump commands up to a specific time, or a statistic computed from recent commands; exercise process data at the specific time; and meal data at the specific time. The one or more observer modules may be configured to receive one or more of the following types of information: a single, processed glucose sample at a specific time, or a history of glucose samples up to a specific time, or a statistic computed from glucose samples up to a specific time; a most recent actual insulin pump command, or a history of recent insulin pump commands up to a specific time, or a statistic computed from recent commands; exercise process data at the specific time; and meal data at the specific time. The one or more observer modules may be configured to process one or more of the following: metabolic measurements; metabolic disturbances; and metabolic treatments.

The metabolic measurements may include one or more of the following: continuous glucose measurements; and insulin measurements. The metabolic disturbances may include one or more of the following: meals; and exercise. The metabolic treatments may include one or more of the following: insulin injections; other pharmaceuticals (hormones) associated with the management of diabetes; treatments for hypoglycemia; and glucagon injections. The other pharmaceuticals may include hormones. The treatments for hypoglycemia may include administering rescue carbohydrates and/or glucagon injections.

The one or more observer modules may be configured to construct and update an internal representation or estimate of a metabolic state of an individual. The one or more observer modules may be configured to transmit metabolic state information to the one or more control modules. The one or more observer modules may be configured to keep an internal representation of a behavioral pattern of an individual. The one or more observer modules may be configured to assess risks of undesirable events. The one or more observer modules may include a short term observer module. The short term observer module may contain information relating to: metabolic state; meal excursion; and/or metabolic state and meal excursion. The short term observer module observes X times per hour, where X is 0<X<7200. It should be appreciated that the frequency may be greater or less as desired or required. The short term observer module may be configured to output one of the following: a vector of estimates of key metabolic states of an individual; a single, processed glucose sample at a specific time, or a history of glucose samples up to a specific time, or a statistic computed from glucose samples up to a specific time; or both the vector of estimates of key metabolic states and the sample or history of processed glucose samples. The vector of estimates may include plasma glucose and plasma insulin.

The one or more observer modules may include a long term observer module. The long term observer module may contain information relating to behavioral profiles. The behavioral profiles may be daily behavioral profiles, but may vary as desired or required. The long term observer module may assess behavioral profiles X times per month, where X is 0<X<60. It should be appreciated that the frequency may be greater or less as desired or required. The long term observer module may be configured to output one or more of the following types of information: a daily profile of an individual's meal behavior as a function of the time of day; a daily profile of the individual's exercise behavior as a function of the time of day; and a daily profile of the individual's utilization of insulin as a function of the time of day. The daily profile of the individual's meal behavior as a function of the time of day may include probabilities of eating at various times throughout the day and/or probabilities of taking various meals defined by carbohydrate, protein, and fat content. The daily profile of the individual's exercise behavior as a function of the time of day may include probabilities of various levels of physical activity throughout the day. The daily profile of the individual's utilization of insulin as a function of the time of day may include a trend analysis for meal and correction boluses and/or a trend analysis for basal rate profiles. The daily profile of an individual's meal behavior may define a probabilistic description of the individual's eating behavior within given meal regimes throughout a day. The behavioral data may be used to define a patient's breakfast regime and then assessing the probability of taking breakfast at any specific point of time within that regime.

The one or more routing modules may include a system coordinator. The system coordinator may be configured to coordinate the distribution of input signals to control one or more modules and allocate different segments of diabetes management to different control modules. The system coordinator may be configured to receive one or more of the following types of information: a single, processed glucose sample at a specific time, or a history of glucose samples up to a specific time, or a statistic computed from glucose samples up to a specific time; a most recent actual insulin pump command, or a history of recent insulin pump commands up to a specific time, or a statistic computed from recent commands; exercise process data at the specific time; meal data at the specific time; a vector of estimates of key metabolic states of an individual; a daily profile of an individual's meal behavior as a function of the time of day; a daily profile of the individual's exercise behavior as a function of the time of day; a daily profile of the individual's utilization of insulin as a function of the time of day; insulin allowed at the specific time by the one or more safety modules as part of a response to meals; insulin allowed at the specific time by the one or more safety modules as part of a response to non-meal disturbances; and estimated glucose excursion due to meals. The system coordinator may be configured to output one or more of the following types of information: a most recent actual insulin pump command; exercise process data at the specific time; meal data at the specific time; a vector of estimates of key metabolic states of an individual; a daily profile of an individual's meal behavior as a function of the time of day; a daily profile of the individual's exercise behavior as a function of the time of day; a daily profile of the individual's utilization of insulin as a function of the time of day; insulin allowed at the specific time by the one or more safety modules as part of a response to meals; insulin allowed at the specific time by the one or more safety modules as part of a response to non-meal disturbances; estimated glucose excursion due to meals; and an estimate of blood glucose concentration offset by the estimated contribution of meals.

The one or more control modules may be configured to perform: daily profile control; compensation control to target; and/or meal control. The daily profile control may include determining a basal rate baseline setting throughout the day, or some other period as desired or required. The compensation control to target may include small adjustments to the basal rate baseline to correct hyperglycemia and reduce the likelihood of hypoglycemia. The meal control may include a schedule of meal insulin following acknowledgement of a particular meal.

The one or more control modules may be configured to include: a first control module to calculate and suggest a basal insulin delivery; a second control module to calculate and suggest compensation of the basal delivery in case of non-meal-related deviations; and a third control module to calculate and suggest meal insulin boluses. The first control module may be configured to receive one or more of the following types of information: a vector of estimates of key metabolic states of an individual; a daily profile of an individual's meal behavior as a function of the time of day; a daily profile of the individual's exercise behavior as a function of the time of day; and a daily profile of the individual's utilization of insulin as a function of the time of day. The first control module may be configured to output a reference insulin infused at a specific time. The second control module may be configured to receive one or more of the following types of information: an estimate of blood glucose concentration offset by the estimated contribution of meals; insulin allowed at the specific time by the one or more safety modules as part of a response to non-meal disturbances; exercise process data at the specific time; a daily profile of the individual's exercise behavior as a function of the time of day; and a daily profile of the individual's utilization of insulin as a function of the time of day. The second control module may be configured to output a recommended residual insulin infusion at a specific time to compensate for non-meal-related disturbances. The third control module may be configured to receive one or more of the following types of information: meal data at the specific time; insulin allowed at the specific time by the one or more safety modules as part of a response to meals; and a daily profile of an individual's meal behavior as a function of the time of day. The third control module may be configured to output one or more of the following types of information: a recommended insulin infusion at a specific time for accommodating meals; and estimated glucose excursion due to meals. The second control module may compensate for a basal delivery up or down. The third control module may be configured to further calculate and suggests pre-meal priming boluses. Moreover, the first control module may operate approximately daily. It should be appreciated that the frequency may be greater or less as desired or required. The second control modules may operate approximately every fifteen to thirty minutes. It should be appreciated that the frequency may be greater or less as desired or required. The third control module may operate approximately every several hours. It should be appreciated that the frequency may be greater or less as desired or required.

The one or more safety modules may be configured to include a system supervision module. The system supervision module may be configured to receive one or more of the following types of information: exercise process data at the specific time; meal data at the specific time; a vector of estimates of key metabolic states of an individual; reference insulin infused at a specific time; a recommended insulin infusion at a specific time for accommodating meals; and a recommended residual insulin infusion at a specific time to compensate for non-meal-related disturbances. The system supervision module may be configured to output one or more of the following types of information: insulin allowed at a specific time by the one or more safety modules as part of a response to meals; insulin allowed at the specific time by the one or more safety modules as part of a response to non-meal disturbances; and total insulin infusion allowed by the safety supervision module at the specific time. The safety supervision module may be configured to: receive information from the one or more observer modules and one or more control modules; determine whether there is an increased risk of hypoglycemia or hyperglycemia; if it is determined that there is an increased risk of hypoglycemia, automatically reduce or discontinue a suggested infusion; and if it is determined that there is an increased risk of hyperglycemia, automatically notify the user of the risk. The increased risk of hypoglycemia may be defined as an increased risk of upcoming hypoglycemia or an increased risk of prolonged hypoglycemia. The suggested infusion may be an insulin infusion or a glucagon infusion. The safety supervision module may further comprise a display device for displaying real-time information to an individual. The individual may be a subject or a doctor, or other user as desired or required. The real-time information may include one or more of the following: a warning of hypoglycemia; a warning of hyperglycemia; a suggestion to reduce insulin delivery; and a suggestion to reject or reduce pre-meal or correction boluses.

Unless defined otherwise, all technical terms used herein have the same meanings as commonly understood by one of ordinary skill in the art of treating diabetes. Specific methods, devices, and materials are described in this application, but any methods and materials similar or equivalent to those described herein can be used in the practice of the present invention. While embodiments of the invention have been described in some detail and by way of exemplary illustrations, such illustration is for purposes of clarity of understanding only, and is not intended to be limiting. Various terms have been used in the description to convey an understanding of the invention; it will be understood that the meaning of these various terms extends to common linguistic or grammatical variations or forms thereof. It will also be understood that when terminology referring to devices, equipment, or drugs has used trade names, brand names, or common names, that these names are provided as contemporary examples, and the invention is not limited by such literal scope. Terminology that is introduced at a later date that may be reasonably understood as a derivative of a contemporary term or designating of a subset of objects embraced by a contemporary term will be understood as having been described by the now contemporary terminology. Further, while some theoretical considerations have been advanced in furtherance of providing an understanding, for example, of the quantitative interrelationships among carbohydrate consumption, glucose levels, and insulin levels, the claims to the invention are not bound by such theory. Moreover, any one or more features of any embodiment of the invention can be combined with any one or more other features of any other embodiment of the invention, without departing from the scope of the invention. Still further, it should be understood that the invention is not limited to the embodiments that have been set forth for purposes of exemplification, but is to be defined only by a fair reading of claims that are appended to the patent application, including the full range of equivalency to which each element thereof is entitled.

Unless clearly specified to the contrary, there is no requirement for any particular described or illustrated activity or element, any particular sequence or such activities, any particular size, speed, material, duration, contour, dimension or frequency, or any particularly interrelationship of such elements. Moreover, any activity can be repeated, any activity can be performed by multiple entities, and/or any element can be duplicated. Further, any activity or element can be excluded, the sequence of activities can vary, and/or the interrelationship of elements can vary. It should be appreciated that aspects of the present invention may have a variety of sizes, contours, shapes, compositions and materials as desired or required.

In summary, while the present invention has been described with respect to specific embodiments, many modifications, variations, alterations, substitutions, and equivalents will be apparent to those skilled in the art. The present invention is not to be limited in scope by the specific embodiment described herein. Indeed, various modifications of the present invention, in addition to those described herein, will be apparent to those of skill in the art from the foregoing description and accompanying drawings. Accordingly, the invention is to be considered as limited only by the spirit and scope of the following claims, including all modifications and equivalents.

Still other embodiments will become readily apparent to those skilled in this art from reading the above-recited detailed description and drawings of certain exemplary embodiments. It should be understood that numerous variations, modifications, and additional embodiments are possible, and accordingly, all such variations, modifications, and embodiments are to be regarded as being within the spirit and scope of this application. For example, regardless of the content of any portion (e.g., title, field, background, summary, abstract, drawing figure, etc.) of this application, unless clearly specified to the contrary, there is no requirement for the inclusion in any claim herein or of any application claiming priority hereto of any particular described or illustrated activity or element, any particular sequence of such activities, or any particular interrelationship of such elements. Moreover, any activity can be repeated, any activity can be performed by multiple entities, and/or any element can be duplicated. Further, any activity or element can be excluded, the sequence of activities can vary, and/or the interrelationship of elements can vary. Unless clearly specified to the contrary, there is no requirement for any particular described or illustrated activity or element, any particular sequence or such activities, any particular size, speed, material, dimension or frequency, or any particularly interrelationship of such elements. Accordingly, the descriptions and drawings are to be regarded as illustrative in nature, and not as restrictive. Moreover, when any number or range is described herein, unless clearly stated otherwise, that number or range is approximate. When any range is described herein, unless clearly stated otherwise, that range includes all values therein and all sub ranges therein. Any information in any material (e.g., a United States/foreign patent, United States/foreign patent application, book, article, etc.) that has been incorporated by reference herein, is only incorporated by reference to the extent that no conflict exists between such information and the other statements and drawings set forth herein. In the event of such conflict, including a conflict that would render invalid any claim herein or seeking priority hereto, then any such conflicting information in such incorporated by reference material is specifically not incorporated by reference herein.

RECITED PUBLICATIONS

The following patents, applications and publications as listed below and throughout this document are hereby incorporated by reference in their entirety herein.

-   1. Bellazzi R, Nucci G, Cobelli C: The subcutaneous route to     insulin-dependent diabetes therapy: closed-loop and partially     closed-loop control strategies for insulin delivery and measuring     glucose concentration. IEEE Eng Med Biol 20: 54-64, 2001. -   2. Bergman R N, Ider Y Z, Bowden C R, Cobelli C. Quantitative     estimation of insulin sensitivity. Am J Physiol. 236:E667-E677,     1979. -   3. Breton M, Kovatchev B P. Analysis, modeling, and simulation of     the accuracy of continuous glucose sensors. J Diabetes Sci Technol     2: 853-862, 2008. -   4. Breton M D, Kovatchev B P, Facchinetti A, Patek S D, Cobelli C.     Adaptive Algorithm Predicting Hypoglycemia in Continuous Glucose     Monitoring (CGM) Data. Proc. 7^(th) Diabetes Technology Meeting, San     Francisco, C A, 2007. -   5. Breton M D, Shields D P, and Kovatchev B P (2008). Optimum     Subcutaneous Glucose Sampling and Fourier Analysis of Continuous     Glucose Monitors. J Diabetes Sci Technol, 2: 495-500. -   6. Buckingham B, Corby E, Clinton P, et al. Preventing Hypoglycemia     Using Predictive Alarm Algorithms and Insulin Pump Suspension,     Diabetes Technol Ther, 11:93-97, 2009. -   7. Chan A, Breton M D and Kovatchev B P. Effects of pulsatile     subcutaneous injections of insulin lispro on plasma insulin     concentration levels. J Diabetes Sci Technol, 2: 844-852, 2008 -   8. Clarke W L and Kovatchev B P. The Artificial Pancreas: How Close     We Are to Closing the Loop? Ped Endocrinol Rev, 4: 314-316, 2007 -   9. Clarke W L and Kovatchev B P. Statistical Tools to Analyze     Continuous Glucose Monitoring Data, Diabetes Technol Ther (in     press—Appendix 3). -   10. Clemens A H, Chang P H, Myers R W. The development of Biostator,     a glucose-controlled insulin infusion system. Horm Metab Res     Supplement, 7: 23-33, 1977. -   11. Cobelli C and Kovatchev B P. Clinical trial of model-predictive     control powered by in silico studies. Proc. 2^(nd) Advanced     Technology & Treatment for Diabetes, Athens, Greece, 2009. -   12. Cobelli C, Ruggeri A: Evaluation of portal/peripheral route and     of algorithms for insulin delivery in the closed-loop control of     glucose in diabetes. A modeling study. IEEE Trans Biomed Eng 30:     93-103, 1983. -   13. Cryer P E. Hypoglycaemia: The limiting factor in the glycaemic     management of type I and type II diabetes. Diabetologia 45: 937-948,     2002. -   14. Dalla Man C, Breton M D, Cobelli C. Physical Activity into the     Meal Glucose-Insulin Model of Type 1 Diabetes: In Silico Studies. J     Diabetes Sci Technol, 3:56-67, 2009. -   15. Dalla Man C, Raimondo D M, Rizza R A, and Cobelli C (2007). GIM,     Simulation Software of Meal Glucose-Insulin Model. J Diabetes     Science and Technology, 1: 323-330. -   16. Dalla Man C, Rizza R A, and Cobelli C (2007). Meal simulation     model of the glucose-insulin system. IEEE Trans Biomed Eng,     54:1740-1749. -   17. Dassau E, Palerm C C, Zisser H, Buckingham B A, Jovanovic L and     F J Doyle III. In silico benchmark platform for artificial     pancreatic b-cell development—a dynamic simulator for closed-loop     control with hardware-in-the-loop. Diabetes Technol Ther 11, 2009 -   18. Dassau E, Zisser H, Palerm C C, Buckingham B A, Jovanovic L,     Doyc III F J. Modular Artificial 1-Cell System: A Prototype for     Clinical Research J Diabetes Sci Technol. 2008; 2:863-872. -   19. Dassau E., Cameron F., Lee H., Bequette B. W., Doyle III F. J.,     GGnter N., Chase H. P. and Buckingham B. A. (2008). “Real-time     Hypoglycemia Prediction Using Continuous Glucose Monitoring (CGM), A     Safety Net to the Artificial Pancreas”. Diabetes 57 (Supl 1): A13.     Proc of the 68th American Diabetes Association Meeting, San     Francisco C A, USA. -   20. Dua P, Doyle F J 3rd, Pistikopoulos E N. “Model-Based Blood     Glucose Control for Type 1 Diabetes via Parametric Programming.”     IEEE Trans Biomed Eng. 53:1478-1491, 2006. -   21. Hovorka R, Canonico V, Chassin L J, et al. Nonlinear model     predictive control of glucose concentration in subjects with type 1     diabetes. Physiol Meas 25:905-920, 2004. -   22. Hovorka R. The future of continuous glucose monitoring: closed     loop. Curr Diabetes Rev 4: 269-79, 2008 -   23. Hovorka R: Continuous glucose monitoring and closed-loop     systems. Diabet Med 23:1-12, 2006. -   24. King C R, Anderson S M, Breton M D, Clarke W L, Kovatchev B P:     Modeling of calibration effectiveness and blood-to-interstitial     glucose dynamics as potential confounders of the accuracy of     continuous glucose sensors during hyperinsulinemic clamp. J Diabetes     Sci Technol, 2007, 1:317-322. -   25. Klonoff D C: The Artificial Pancreas: How Sweet Engineering Will     Solve Bitter Problems. J Diabetes Sci Technol, 1: 72-81, 2007. -   26. Kovatchev B P and Clarke W L. Peculiarities of the Continuous     Glucose Monitoring Data Stream and Their Impact on Developing     Closed-Loop Control Technology. J Diabetes Sci Technol, 2:158-163,     2008. -   27. Kovatchev B P et al. Personalized subcutaneous model-predictive     closed-loop control of T1DM: Pilot studies in the USA and Italy,     Proc 69^(th) ADA Scientific Session, 228-OR (in press—Appendix 2). -   28. Kovatchev B P et al. Study of Closed-Loop Glucose Control in     Type 1 Diabetes. Food and Drug Administration Investigational Device     Exemption G080048, 2008. -   29. Kovatchev B P, Anderson S M, Heinemann L, Clarke W L. Comparison     of the numerical and clinical accuracy of four continuous glucose     monitors. Diabetes Care, 31: 1160-1164, 2008 -   30. Kovatchev B P, Anderson S M, Otto E. Field Glucose Variability     Index is Related to Laboratory Measures of Insulin Sensitivity and     Hypoglycemia Counterregulation. Diabetologia, 49, Supplement 1:     A855, 2005. -   31. Kovatchev B P, Breton M D, Dalla Man C, Cobelli C. In Silico     Preclinical Trials: A Proof of Concept in Closed-Loop Control of     Type 1 Diabetes. J Diabetes Sci Technol 3: 44-55, 2009. -   32. Kovatchev B P, Breton M D, Dalla Man C, Cobelli C. In Silico     model and computer simulation environment approximating the human     glucose/insulin utilization. Food and Drug Administration Master     File MAF 1521, 2008. -   33. Kovatchev B P, Breton M D, Dalla Man C, Cobelli C. In Silico     Preclinical Trials: A Proof of Concept in Closed-Loop Control of     Type 1 Diabetes. J Diabetes Sci Technol, 2009, 3:44-55. -   34. Kovatchev B P, Clarke W L, Breton M, Brayman K, McCall A.     Quantifying Temporal Glucose Variability in Diabetes via Continuous     Glucose Monitoring: Mathematical Methods and Clinical Application.     Diabetes Tech Ther, 7: 849-862, 2005. -   35. Kovatchev B P, Cox D J, Farhy L S, Straume M, Gonder-Frederick L     A, Clarke, W L. Episodes of Severe Hypoglycemia in Type 1 Diabetes     are Preceded, and Followed, within 48 Hours by Measurable     Disturbances in Blood Glucose. J of Clinical Endocrinology and     Metabolism, 85: 4287-4292, 2000. -   36. Kovatchev B P, Cox D J, Gonder-Frederick L A and W L Clarke.     Symmetrization of the blood glucose measurement scale and its     applications. Diabetes Care 20: 1655-1658, 1997 -   37. Kovatchev B P, Cox D J, Gonder-Frederick L A Young-Hyman D,     Schlundt D, Clarke W L. Assessment of risk for severe hypoglycemia     among adults with IDDM: Validation of the Low Blood Glucose Index,     Diabetes Care 21: 1870-1875, 1998. -   38. Kovatchev B P, Cox D J, Kumar A, Gonder-Frederick L A and W L     Clarke. Algorithmic Evaluation of Metabolic Control and Risk of     Severe Hypoglycemia in Type 1 and Type 2 Diabetes Using     Self-Monitoring Blood Glucose (SMBG) Data. Diabetes Technol Ther, 5:     817-828, 2003 -   39. Kovatchev B P, Cox D J, Straume M, Farhy L S. Association of     Self-monitoring Blood Glucose Profiles with Glycosylated Hemoglobin     (2000). In: Methods in Enzymology, vol. 321: Numerical Computer     Methods, Part C: 410-417 M. Johnson & L. Brand, Eds., Academic     Press, N Y. -   40. Kovatchev B P, Gonder-Frederick L A, Cox D J, Clarke W L:     Evaluating the accuracy of continuous glucose monitoring sensors:     Continuous Glucose Error-Grid Analysis (C G-EGA) illustrated by     Therasense Freestyle Navigator™ data. Diabetes Care 27:1922-1928,     2004. -   41. Kovatchev B P, Hovorka R, Weinzimer S for the JDRF Artificial     Pancreas Consortium. Outcome measures: Judging the effectiveness of     closed-loop control. JDRF AP Consortium Position Paper. -   42. Kovatchev B P, Otto E, Cox D J, Gonder-Frederick L A, Clarke     W L. Evaluation of a New Measure of Blood Glucose Variability in     Diabetes. Diabetes Care, 29: 2433-2438, 2006. -   43. Kovatchev B P, Straume M, Cox D J, Farhy L S. Risk analysis of     blood glucose data: A quantitative approach to optimizing the     control of insulin dependent diabetes. J of Theoretical Medicine,     3:1-10, 2001. -   44. Magni L, Raimondo F, Bossi L, Dalla Man C, De Nicolao G,     Kovatchev B P, Cobelli C. Model Predictive Control of Type 1     Diabetes: An In Silico Trial J Diabetes Sci Technol, 1:804-812, 2007 -   45. Magni L, Raimondo F, Dalla Man C, Breton M D, Patek S, De     Nicolao G, Cobelli C, and Kovatchev B P. Evaluating the Efficacy of     Closed-Loop Glucose Regulation via Control-Variability Grid     Analysis. J Diabetes Sci Technol, 2: 630-635, 2008 -   46. McCall A, Cox D J, Crean J, Gloster M, and Kovatchev B P. A     Novel Analytical Method for Assessing Glucose Variability: Using     CGMS in Type 1 Diabetes Mellitus. Diabetes Technology and     Therapeutics, 8: 644-653, 2006. -   47. Nucci G., Cobelli C. Models of subcutaneous insulin kinetics. A     critical review. Comput Methods Programs Biomed., 62:249-57 Review,     2000. -   48. Owens C, Zisser H, Jovanovic L, Srinivasan B, Bonvin D, Doyle F     J 3rd. Run-to-run control of blood glucose concentrations for people     with Type 1 diabetes mellitus. IEEE Trans Biomed Eng. 53:996-1005,     2006 -   49. Palerm C C, Zisser H, Bevier W C, Jovanovic L, Doyle F J 3rd.     Prandial insulin dosing using run-to-run control: application of     clinical data and medical expertise to define a suitable performance     metric. Diabetes Care. 30:1131-1136, 2007 -   50. Parker R S, Doyle F J 3rd, Peppas N A. A model-based algorithm     for blood glucose control in Type I diabetic patients. IEEE Trans     Biomed Eng, 48:148-157, 1999. -   51. Patek S D, Breton M D, Chen Y, Solomon C, and Kovatchev B P.     LQG-Based Closed-Loop Control of Type 1 Diabetes. J Diabetes Sci     Technol, 1: 834-841, 2007. -   52. Patek S D, Breton M D, Cobelli C, Dalla Man C, and Kovatchev     B P. Adaptive Meal Detection Algorithm Enabling Closed-Loop Control     in Type 1 Diabetes. Proc. 7^(th) Diabetes Technology Meeting, San     Francisco, CA, 2007. -   53. Patek S D, Breton M D, Hughes C, and Kovatchev B P. Control of     Hypoglycemia via Estimation of Active Insulin, Glucose Forecasts,     and Risk-Based Insulin Reduction. Proc. 2^(nd) Advanced Technology &     Treatment for Diabetes, Athens, Greece, 2009. -   54. Pillonetto G, Caumo A, Sparacino G, Cobelli C. A new dynamic     index of insulin sensitivity. IEEE Transactions on Biomedical     Engineering. 53: 369-79, 2006. -   55. Sorensen J T: A Physiologic Model of Glucose Metabolism in Man     and its Use to Design and Assess Improved Insulin Therapies for     Diabetes, Ph.D. dissertation, Dept Chemical Engineering, MIT, 1985. -   56. Sparacino G, Zanderigo F, Corazza G, Maran A, Facchinetti A,     Cobelli C. Glucose concentration can be predicted ahead in time from     continuous glucose monitoring sensor time-series. IEEE Trans. Biomed     Eng, 54:931-937, 2007 -   57. Steil G M, Rebrin K, Darwin C, Hariri F, Saad M F. Feasibility     of automating insulin delivery for the treatment of type 1 diabetes.     Diabetes 55: 3344-3350, 2006. -   58. The Diabetes Control and Complications Trial Research Group. The     effect of intensive treatment of diabetes on the development and     progression of long-term complications of insulin-dependent diabetes     mellitus. N Engl. J Med 329: 978-986, 1993. -   59. U K Prospective Diabetes Study Group (UKPDS). Intensive     blood-glucose control with sulphonylureas or insulin compared with     conventional treatment and risk of complications in patients with     type 2 diabetes. Lancet 352: 837-853, 1998. -   60. Weinzimer S A, Steil G M, Swan K L, Dziura J, Kurtz N,     Tamborlane W V: Fully automated closed-loop insulin delivery versus     semi-automated hybrid control in pediatric patients with type 1     diabetes using an artificial pancreas. Diabetes Care 31:934-939,     2008. -   61. Zanderigo F, Sparacino G, Kovatchev B, Cobelli C. Glucose     Prediction Algorithms from continuous monitoring: Assessment of     accuracy via Continuous Glucose-Error Grid Analysis. J Diabetes Sci     Technol, 1: 645-651, 2007 -   62. Zisser H, Jovanovic L, Doyle F J 3rd, Ospina P, Owens C.     Run-to-Run Control of Meal-Related Insulin Dosing. Diab Technol     Ther. 7:48-57, 2005

REFERENCES

The devices, structures, systems, computer program products, and methods of various embodiments of the invention disclosed herein may utilize aspects disclosed in the following references, applications, publications and patents and which are hereby incorporated by reference herein in their entirety:

-   -   A. International Patent Application Serial No.         PCT/US2010/025405, entitled “Method, System and Computer Program         Product for CGM-Based Prevention of Hypoglycemia via         Hypoglycemia Risk Assessment and Smooth Reduction Insulin         Delivery,” filed Feb. 25, 2010     -   B. International Patent Application Serial No.         PCT/US2009/065725, filed Nov. 24, 2009, entitled “Method,         System, and Computer Program Product for Tracking of Blood         Glucose Variability in Diabetes from Data,”     -   C. PCT/US2008/082063, entitled “Model Predictive Control Based         Method for Closed-Loop Control of Insulin Delivery in Diabetes         Using Continuous Glucose Sensing”, filed Oct. 31, 2008.     -   D. PCT/US2008/069416, entitled “Method, System and Computer         Program Product for Evaluation of Insulin Sensitivity,         Insulin/Carbohydrate Ratio, and Insulin Correction Factors in         Diabetes from Self-Monitoring Data”, filed Jul. 8, 2008.     -   E. PCT/US2008/067725, entitled “Method, System and Computer         Simulation Environment for Testing of Monitoring and Control         Strategies in Diabetes,” filed Jun. 20, 2008.     -   F. PCT/US2008/067723, entitled “LQG Artificial Pancreas Control         System and Related Method”, filed on Jun. 20, 2008.     -   G. U.S. Ser. No. 12/516,044, filed May 22, 2009, entitled         “Method, System, and Computer Program Product for the Detection         of Physical Activity by Changes in Heart Rate, Assessment of         Fast Changing Metabolic States, and Applications of Closed and         Open Control Loop in Diabetes;”     -   H. PCT/US2007/085588 not yet published filed Nov. 27, 2007,         entitled “Method, System, and Computer Program Product for the         Detection of Physical Activity by Changes in Heart Rate,         Assessment of Fast Changing Metabolic States, and Applications         of Closed and Open Control Loop in Diabetes;”     -   I. U.S. Ser. No. 11/943,226, filed Nov. 20, 2007, entitled         “Systems, Methods and Computer Program Codes for Recognition of         Patterns of Hyperglycemia and Hypoglycemia, Increased Glucose         Variability, and Ineffective Self-Monitoring in Diabetes;”     -   J. U.S. patent application Ser. No. 11/578,831, filed Oct. 18,         2006 entitled “Method, System and Computer Program Product for         Evaluating the Accuracy of Blood Glucose Monitoring         Sensors/Devices”.     -   K. PCT International Application Serial No. PCT/US2005/013792,         filed Apr. 21, 2005, entitled “Method, System, and Computer         Program Product for Evaluation of the Accuracy of Blood Glucose         Monitoring Sensors/Devices;”     -   L. PCT International Application Serial No. PCT/US01/09884,         filed Mar. 29, 2001, entitled “Method, System, and Computer         Program Product for Evaluation of Glycemic Control in Diabetes         Self-Monitoring Data;”     -   M. U.S. Pat. No. 7,025,425 B2 issued Apr. 11, 2006, entitled         “Method, System, and Computer Program Product for the Evaluation         of Glycemic Control in Diabetes from Self-Monitoring Data;”     -   N. U.S. patent application Ser. No. 11/305,946 filed Dec. 19,         2005 entitled “Method, System, and Computer Program Product for         the Evaluation of Glycemic Control in Diabetes from         Self-Monitoring Data” (Publication No. 2006/0094947);     -   O. PCT International Application Serial No. PCT/US2003/025053,         filed Aug. 8, 2003, entitled “Method, System, and Computer         Program Product for the Processing of Self-Monitoring Blood         Glucose (SMBG) Data to Enhance Diabetic Self-Management;”     -   P. U.S. patent application Ser. No. 10/524,094 filed Feb. 9,         2005 entitled “Managing and Processing Self-Monitoring Blood         Glucose” (Publication No. 2005/214892);     -   Q. U.S. Ser. No. 12/065,257, filed Aug. 29, 2008, entitled         “Accuracy of Continuous Glucose Sensors;”     -   R. PCT International Application Serial No PCT/US2006/033724,         filed Aug. 29, 2006, entitled “Method for Improvising Accuracy         of Continuous Glucose Sensors and a Continuous Glucose Sensor         Using the Same;”     -   S. U.S. Ser. No. 12/159,891, filed Jul. 2, 2008, entitled         “Method, System and Computer Program Product for Evaluation of         Blood Glucose Variability in Diabetes from Self-Monitoring         Data;”     -   T. PCT International Application No. PCT/US2007/000370, filed         Jan. 5, 2007, entitled “Method, System and Computer Program         Product for Evaluation of Blood Glucose Variability in Diabetes         from Self-Monitoring Data;”     -   U. U.S. patent application Ser. No. 11/925,689 and PCT         International Patent Application No. PCT/US2007/082744, both         filed Oct. 26, 2007, entitled “For Method, System and Computer         Program Product for Real-Time Detection of Sensitivity Decline         in Analyte Sensors;”     -   V. U.S. Ser. No. 10/069,674, filed Feb. 22, 2002, entitled         “Method and Apparatus for Predicting the Risk of Hypoglycemia;”     -   W. PCT International Application No. PCT/US00/22886, filed Aug.         21, 2000, entitled “Method and Apparatus for Predicting the Risk         of Hypoglycemia;”     -   X. U.S. Pat. No. 6,923,763 B1, issued Aug. 2, 2005, entitled         “Method and Apparatus for Predicting the Risk of Hypoglycemia;”     -   Y. U.S. Patent Application No. US 2004/0254434 A1, “Glucose         Measuring Module and “Insulin Pump Combination”, Dec. 16, 2004.     -   Z. U.S. Patent Application Publication No. US 2009/00697456 A1,         Estes, et al., “Operating an Infusion Pump System”, Mar. 12,         2009. 

1-9. (canceled)
 10. A method of adaptive diabetes control by a diabetes control system, the method comprising: receiving a first glucose input reflective of glucose sample processing by a continuous glucose monitor (CGM) sensor; receiving a first insulin delivery command, the first insulin delivery command representing a most recent command to an insulin pump relative to a time of the first glucose input; computing, in a first control function, a first reference insulin infusion based on the first glucose input, the first insulin delivery command, and a first profile descriptive of a statistical insulin utilization; computing, in a second control function, a first residual insulin infusion based on the first glucose input and a second profile descriptive that is at least a function of time of day; computing, in a third control function, a first meal-accommodation insulin infusion based on a third profile of a meal behavior; and determining a first total insulin infusion based on the first reference insulin infusion, the first residual insulin infusion, and the first meal-accommodation insulin infusion.
 11. The method of claim 10, comprising the diabetes control system: commanding the insulin pump to deliver the determined first total insulin infusion via a second insulin delivery command.
 12. The method of claim 11, wherein the diabetes control system is external to the insulin pump and the commanding comprises: sending the second insulin delivery command to the insulin pump over a wireless communications channel.
 13. The method of claim 11, wherein the computing a first residual insulin infusion comprises: managing, in the second control function, a blood glucose excursion not accounted for by the third control function.
 14. The method of claim 11, comprising, by the diabetes control system: estimating, in the third control function, a first glucose excursion from a target due to at least one of a meal or a response to a meal; and estimating, based on the estimated first glucose excursion, a second glucose excursion from the target that is not due to the at least one of a meal or a response to a meal, wherein the first residual insulin infusion is computed based on the estimated second glucose excursion.
 15. The method of claim 11, wherein the first reference insulin infusion comprises: a suggested basal delivery, the first residual insulin infusion compensating the suggested basal delivery either upward to correct hyperglycemia or downward to reduce a likelihood of hypoglycemia.
 16. The method of claim 11, the method comprising, by the diabetes control system: receiving a second glucose input reflective of glucose sample processing by the CGM sensor; receiving the second insulin delivery command, the received second insulin delivery command representing a most recent command to the insulin pump relative to a time of the second glucose input; computing, in the first control function, a second reference insulin infusion based on the second glucose input, the received second insulin delivery command, and the first profile descriptive of the statistical insulin utilization; computing, in the second control function, a second residual insulin infusion based on the second glucose input and the second profile descriptive that is at least a function of time of day; computing, in the third control function, a second meal-accommodation insulin infusion based on the third profile of meal behavior; determining a second total insulin infusion based on the second reference insulin infusion, the second residual insulin infusion, and the second meal-accommodation insulin infusion; and commanding the insulin pump to deliver the determined second total insulin via a third insulin delivery command.
 17. The method of claim 11, comprising the diabetes control system: assessing at least one of the first profile, the second profile, or the third profile and, subsequent to the assessing: receiving a second glucose input reflective of glucose sample processing by the CGM sensor; receiving the second insulin delivery command, the received second insulin delivery command representing a most recent command to the insulin pump relative to a time of the second glucose input; computing, in the first control function, a second reference insulin infusion based on the second glucose input, the received second insulin delivery command, and the first profile descriptive statistical insulin utilization; computing, in the second control function, a second residual insulin infusion based on the second glucose input and the second profile descriptive that is at least a function of time of day; computing, in the third control function, a second meal-accommodation insulin infusion based on the third profile of the meal behavior; determining a second total insulin infusion based on the second reference insulin infusion, the second residual insulin infusion, and the second meal-accommodation insulin infusion; and commanding the insulin pump to deliver the determined second total insulin infusion via a third insulin delivery command.
 18. The method of claim 11, wherein the first residual insulin infusion is computed via at least one of the following: a model predictive control (MPC) technique; a dynamic model of interactions of a glucose target and residual insulin; a linear quadratic Gaussian (LQG) control methodology; a learning model predictive control (LMPC) methodology; and/or a proportional-integral-derivative (PID) control methodology.
 19. The method of claim 10, comprising the diabetes control system: assessing the first profile, the second profile, and the third profile at least monthly.
 20. A system for diabetes control, the system comprising: at least one processor; and memory having control logic stored thereon that is executable by the at least one processor to perform operations including: receiving a first glucose input reflective of glucose sample processing by a continuous glucose monitor (CGM) sensor configured for a patient; receiving a first insulin delivery command, the first insulin delivery command representing a most recent command to an insulin pump relative to a time of the first glucose input; computing, in a first control function, a first reference insulin infusion based on the first glucose input, the first insulin delivery command, and a first profile descriptive of a statistical insulin utilization; computing, in a second control function, a first residual insulin infusion based on the first glucose input and a second profile descriptive that is at least a function of time of day; computing, in a third control function, a first meal-accommodation insulin infusion based on a third profile of meal behavior; and determining a first total insulin infusion based on the first reference insulin infusion, the first residual insulin infusion, and the first meal-accommodation insulin infusion.
 21. The system of claim 20, wherein the operations comprise: commanding the insulin pump to deliver the determined first total insulin infusion via a second insulin delivery command.
 22. The system of claim 21, wherein the computing a first residual insulin infusion comprises: managing, in the second control function, a blood glucose excursion not accounted for by the third control function.
 23. The system of claim 21, wherein the operations comprise: estimating, in the third control function, a first glucose excursion from a target due to at least one of a meal or a response to a meal; and estimating, based on the estimated first glucose excursion, a second glucose excursion from the target that is not due to the at least one of a meal or a response to a meal, wherein the first residual insulin infusion is computed based on the estimated second glucose excursion.
 24. The system of claim 21, wherein the first reference insulin infusion includes a suggested basal delivery, the first residual insulin infusion compensating the suggested basal delivery either upward to correct hyperglycemia or downward to reduce a likelihood of hypoglycemia.
 25. The system of claim 21, wherein the operations comprise: receiving a second glucose input reflective of glucose sample processing by the CGM sensor; receiving the second insulin delivery command, the received second insulin delivery command representing a most recent command to the insulin pump relative to a time of the second glucose input; computing, in the first control function, a second reference insulin infusion based on the second glucose input, the received second insulin delivery command, and the first profile descriptive of the statistical insulin utilization; computing, in the second control function, a second residual insulin infusion based on the second glucose input and the second profile descriptive that is at least a function of time of day; computing, in the third control function, a second meal-accommodation insulin infusion based on the third profile of the meal behavior; determining a second total insulin infusion based on the second reference insulin infusion, the second residual insulin infusion, and the second meal-accommodation insulin infusion; and commanding the insulin pump to deliver the determined second total insulin infusion via a third insulin delivery command.
 26. The system of claim 21, wherein the operations comprise: assessing at least one of the first profile, the second profile, or the third profile and, subsequent to the assessing: receiving a second glucose input reflective of glucose sample processing by the CGM sensor; receiving the second insulin delivery command, the received second insulin delivery command representing a most recent command to the insulin pump relative to a time of the second glucose input; computing, in the first control function, a second reference insulin infusion based on the second glucose input, the received second insulin delivery command, and the first profile descriptive of the statistical insulin utilization; computing, in the second control function, a second residual insulin infusion based on the second glucose input and the second profile descriptive of the patient at least a function of time of day; computing, in the third control function, a second meal-accommodation insulin infusion based on the third profile of the meal behavior; determining a second total insulin infusion based on the second reference insulin infusion, the second residual insulin infusion, and the second meal-accommodation insulin infusion; and commanding the insulin pump to deliver the determined second total insulin infusion via a third insulin delivery command.
 27. The system of claim 21, wherein the operations comprise: assessing the first profile, the second profile, and the third profile at least monthly.
 28. The system of claim 21, wherein the system is a personal device comprising: a display and an interface for wireless communication, and the operations comprise: sending the second insulin delivery command to the insulin pump via the interface for wireless communication.
 29. The system of claim 21, wherein the system is a personal device comprising: a display and an interface for cellular communication.
 30. A computer-program product having a non-transitory computer-usable medium with computer program logic stored thereon, the computer program logic being configured to be executed to perform operations comprising: receiving a first glucose input reflective of glucose sample processing by a continuous glucose monitor (CGM) sensor configured for a patient; receiving a first insulin delivery command, the first insulin delivery command representing a most recent command to an insulin pump relative to a time of the first glucose input; computing, in a first control function, a first reference insulin infusion based on the first glucose input, the first insulin delivery command, and a first profile descriptive of a statistical insulin utilization; computing, in a second control function, a first residual insulin infusion based on the first glucose input and a second profile descriptive that is at least a function of time of day; computing, in a third control function, a first meal-accommodation insulin infusion based on a third profile of meal behavior; and determining a first total insulin infusion based on the first reference insulin infusion, the first residual insulin infusion, and the first meal-accommodation insulin infusion.
 31. The computer-program product of claim 30, wherein the operations comprise: commanding the insulin pump to deliver the determined first total insulin via a second insulin delivery command.
 32. The computer-program product of claim 22, wherein the computing a first residual insulin infusion comprises: managing, in the second control function, a blood glucose excursion not accounted for by the third control function.
 33. The computer-program product of claim 31, wherein the operations comprise: estimating, in the third control function, a first glucose excursion from a target due to at least one of a meal or a response to a meal; and estimating, based on the estimated first glucose excursion, a second glucose excursion from the target that is not due to the at least one of a meal or a response to a meal, wherein the first residual insulin infusion is computed based on the estimated second glucose excursion.
 34. The computer-program product of claim 31, wherein the first reference insulin infusion comprises: a suggested basal delivery, the first residual insulin infusion compensating the suggested basal delivery either upward to correct hyperglycemia or downward to reduce a likelihood of hypoglycemia.
 35. The computer-program product of claim 31, wherein the operations comprise: receiving a second glucose input reflective of glucose sample processing by the CGM sensor; receiving the second insulin delivery command, the received second insulin delivery command representing a most recent command to the insulin pump relative to a time of the second glucose input; computing, in the first control function, a second reference insulin infusion based on the second glucose input, the received second insulin delivery command, and the first profile descriptive of the statistical insulin utilization; computing, in the second control function, a second residual insulin infusion based on the second glucose input and the second profile descriptive that is at least a function of time of day; computing, in the third control function, a second meal-accommodation insulin infusion based on the third profile of the meal behavior; determining a second total insulin infusion based on the second reference insulin infusion, the second residual insulin infusion, and the second meal-accommodation insulin infusion; and commanding the insulin pump to deliver the determined second total insulin infusion via a third insulin delivery command.
 36. The computer-program product of claim 31, wherein the operations comprise: assessing at least one of the first profile, the second profile, or the third profile and, subsequent to the assessing: receiving a second glucose input reflective of glucose sample processing by the CGM sensor; receiving the second insulin delivery command, the received second insulin delivery command representing a most recent command to the insulin pump relative to a time of the second glucose input; computing, in the first control function, a second reference insulin infusion based on the second glucose input, the received second insulin delivery command, and the first profile descriptive of the statistical insulin utilization; computing, in the second control function, a second residual insulin infusion based on the second glucose input and the second profile descriptive that is at least a function of time of day; computing, in the third control function, a second meal-accommodation insulin infusion based on the third profile of the meal behavior; determining a second total insulin infusion based on the second reference insulin infusion, the second residual insulin infusion, and the second meal-accommodation insulin infusion; and commanding the insulin pump to deliver the determined second total insulin infusion to the patient via a third insulin delivery command.
 37. The computer-program product of claim 31, wherein the operations comprise: assessing the first profile, the second profile, and the third profile at least monthly.
 38. The computer-program product of claim 31, wherein the computer program logic is implemented on a personal device having a display and an interface for wireless communication, and the operations comprise: sending the second insulin delivery command to the insulin pump via the interface for wireless communication.
 39. The computer-program product of claim 31, wherein the computer program logic is implemented on a personal device having a display and an interface for cellular communication. 